Constituent Gluon Content of the Static Quark- Antiquark State in Coulomb Gauge 
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Motivated by the gluon-chain model of flux tube formation, we compute and diagonalize the transfer matrix 
in lattice SU(2) gauge theory for states containing heavy static quark-antiquark sources, with separations up to 
one fermi. The elements of the transfer matrix are calculated by variational Monte Carlo methods, in a basis of 
states obtained by acting on the vacuum state with zero, one, and two-gluon operators in Coulomb gauge. The 
color Coulomb potential is obtained from the zero gluon to zero gluon element of the transfer matrix, and it is 
well-known that while this potential is asymptotically linear, it has a slope which is two to three times larger 
than the standard asymptotic string tension. We show that the addition of one and two gluon states results in 
a potential which is still linear, but the disagreement with the standard asymptotic string tension is reduced to 
38% at the largest lattice coupling we have studied. 



I. INTRODUCTION 



It is well known that the color Coulomb potential in non- 
abelian pure gauge theories is asymptotically linear, but that 
the corresponding string tension O cou i is greater than the stan- 
dard string tension by anywhere from 100-200%, depend- 
ing on the lattice coupling and the gauge group [1, 2]. This 
discrepancy should not be too surprising, since the color 
Coulomb potential is only an upper bound on the static quark 
potential [3], but the more relevant fact is that the Coulombic 
force arises from (dressed) one-gluon exchange between static 
sources. A serious criticism of one-gluon exchange models 
of the confining potential, whether in Coulomb or covariant 
gauges, is that they provide no explanation whatever for the 
collimation of color electric fields into flux tubes, and, as a 
consequence, no obvious reason for the absence of long-range 
dipole forces. Since long-range dipole forces are incompatible 
with a massive spectrum, a theory of confinement via dressed 
one-gluon exchange is in danger of inconsistency, as well as 
conflicting with the numerical evidence for color-electric flux 
tubes. 

There is, however, a picture of how flux tubes can form in 
Coulomb gauge, in which strong one-particle exchange forces 
are an important ingredient. This is the "gluon-chain" model, 
put forward by Thorn and one of the present authors [4]. The 
idea is that as a heavy quark and and antiquark separate, they 
pull out a sequence of constituent gluons between them, as il- 
lustrated in Fig. 1. The constituent gluons are bound together 
by Coulombic nearest-neighbor interactions, and the ensem- 
ble of gluons resembles a discretized string. In the original 
proposal it was supposed that the Coulombic force at some 
distance scale rises faster than linearly, and that the role of the 
constituent gluons is to prevent a corresponding faster-than- 
linear rise in the static quark potential, by effectively plac- 
ing an upper limit on color charge separation. We now know 
that the Coulomb potential is itself linear at large scales. So 
the role of constituent gluons must be to somehow reduce the 
magnitude of the static quark potential, without at the same 
time destroying the linearity property. That such an effect 
might take place was shown analytically by Krupinski and 
Szczepaniak [5], in the context of a particular proposal for 



the form of the Coulomb gauge vacuum state. It is desirable, 
however, to see if constituent gluons can have this conjectured 
effect without making any assumptions about the form of the 
vacuum. That goal, at present, can only be achieved numeri- 
cally, via lattice Monte Carlo simulations. 

There has been only one attempt [6], carried out twenty 
years ago, to study the constituent gluon content of the QCD 
flux tube in Coulomb gauge by numerical simulation, and 
thereby test the gluon chain model. Our intention in the 
present article is to greatly improve on that earlier effort, 
and to quantify the reduction in the energy of a static quark- 
antiquark state which can be achieved with a handful of con- 
stituent gluon operators. 
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FIG. 1: The gluon chain model. As a quark-antiquark pair moves 
apart, they pull out a chain of constituent gluons between them. 
Open/filled shadings of the constituent gluons represent the two ma- 
trix indices of the A?* operator. Dashed lines indicate both repeated 
matrix indices, and nearest-neighbor Coulomb interactions. 



II. THE TRANSFER MATRIX IN A VARIATIONAL BASIS 

The Euclidean-time evolution operator in lattice gauge the- 
ory, in some physical gauge such as Coulomb gauge, is the 
transfer matrix 



T = exp[— Ha) 



(2.1) 
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where a is the lattice spacing and H the Hamiltonian. It is 
useful to consider the rescaled operator 



1 



(^olTl^o) 
exp[— (H — Eo)a] 



(2.2) 



where 'I'o, Eq, are the ground state and vacuum energy, re- 
spectively. To compute the static quark potential of a quark- 
antiquarkpair separated by a distance R, one would ideally di- 
agonalize the transfer matrix in the infinite-dimensional sub- 
space of states which contain a single massive quark, and 
a single massive antiquark, located at sites x and y with 
R = |x — yj. The minimal energy eigenstate of the transfer 
matrix, in this subspace, is the state with the largest eigen- 
value X max of T, and the static quark potential, in lattice units, 
is given by 



V(R) 



\og(X inax ) 



(2.3) 



(from here on we will work in lattice units). If we consider 
any two states \a) and \b) which are obtained by acting on the 
vacuum state with operators Q a and Qt,, respectively, then 



(a\T\b) = (Ql(t + l)Q b (t)) 



(2.4) 



where the notation Q(t) indicates that the operator is to be 
evaluated using link variables on a hypersurface of fixed time 
t . For given operators Q a ,b, the rhs of eq. (2.4) can be evalu- 
ated by lattice Monte Carlo simulation. 

If the quark and antiquark are not too far apart, then accord- 
ing to the gluon-chain picture it should be possible to express 
the minimal energy state in terms of the vacuum state ^o, and 
a handful of operators operators acting on that state. These op- 
erators create a massive antiquark at lattice site x, a massive 
quark at site y, and a small number of constituent gluons. Let 
us denote the gluonic operators by Qk, deferring their actual 
construction to the next section, with corresponding states 

|*)=?'(x)GgV(y)| , Fo) k=l,2,-,M (2.5) 

The {Qk} are functionals of the lattice gauge field, and depend 
on x, y, and some number of variational parameters. In our 
construction, described below in section III, we use a single 
variational parameter, denoted p. One can then compute, by 
lattice Monte Carlo simulation, the quantities 

O mn = (m\n) 

= (^v[Ql(t)Q n (t)]) 
t, nn = (m\T\n) 

= (^[Ql(t + i)U f MQn(t)U (x L ,t)}) (2.6) 



We cannot diagonalize the transfer matrix T in the basis 
{\k}}, because these states are not orthogonal, in general, nor 
are they normalized. So the next thing to do is to construct 
a new set of orthogonal states {|&')} from the non-orthogonal 



states by the Gram-Schmidt orthogonalization procedure 

= LD kj \j) (2.7) 



Defining 



Nj = (/|/) 

= EE^w^H") 



(2.8) 



we have 



I* 7 ) = l*>-I^I^«H*>I>»l«} 



j=l N .i m 



jfe-1 
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D: 



E Tf E DP'Mk) 

j=l 



= 1 



In) (2.9) 



The D mn coefficients can be determined iteratively by the fol- 
lowing set of relations: 

-Z k j^Nj l D jn lj m ^D jm (m\k) n<k 

1 n = k (2.10) 

n > k 

Normalize, 



(2.11) 



and calculate the M x M matrix 

Tn = (fo|r|ty) 



D in Dj m (n\T\m) 



1 



E E d*pj 



(2.12) 



Note that this matrix is derived from the quantities O mn and 
t mn of eq. (2.6), both of which are calculated by lattice Monte 
Carlo. 

The final step is to extract the largest eigenvalue X max of 
the finite matrix 7y. Since this is a variational approach, the 
computation has to be repeated for a variety of values of the 
variational parameter p, at each quark-antiquark separation R. 
The value of p which minimizes — log(X max ) gives the best 
variational estimate 

V chain (R) = -\og(\„ wx ) (2.13) 

for the static quark potential. We will refer to this variational 
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result as the "gluon-chain potential", although in practice the 
variational state will be a superposition of states containing at 
most two constituent gluons. 



HI. SMEARED ONE AND TWO-GLUON OPERATORS 

In order to actually carry out the procedure described above, 
we must specify a set of operators {Qi } which, operating on 
the Coulomb gauge vacuum, generate a set of variational trial 
states {\k)}. 

A gluon chain consists of a number of constituent gluons 
lying between the heavy sources, with the ordering of the glu- 
ons in their color-indices correlating with the spatial positions 
of gluons between the static sources. Suppose, e.g., that the 
static quark-antiquark sources are located on the x-axis at lo- 
cations x = and x~R. Constituent gluons which are located 
well outside this interval in their x-coordinates, or which have 
large transverse displacements away from the x-axis, will be 
costly in terms of Coulombic interaction energy (recall that 
this energy itself increases linearly with color charge separa- 
tion). On the other hand, highly localized constituent gluons 
are costly in terms of kinetic energy. As usual, the variational 
procedure attempts to find an optimal compromise between 
kinetic energy, which favors spatial derealization, and inter- 
action energy, which in this case favors small transverse dis- 
placement from the line joining the quark sources. Derealiza- 
tion in the x-direction is achieved by a superposition of gluon 
operators (in the Qk) at different locations along the x 6 [0,R] 
interval. 1 Derealization in the transverse directions can be 
obtained by constructing A-field operators on the lattice, in 
which the high-frequency components of the A-field, in the 
transverse directions, are gaussian suppressed. 

The delocalized, or "transverse-smoothed", operators 
A? (x,t,j) are constructed in the following way: The L 4 lattice 
is fixed to Coulomb gauge by standard methods (simulated 
annealing + over-relaxation), and the usual lattice A-field vari- 
able, in the SU(2) gauge group, is defined in terms of the link 
variables via 



At(x,f) 



±(u k (x,t)-Ul(x,t] 



Denote the finite Fourier transform of this array by 
a(n\,n2,ri3), with indices running over the same range of 
to L — 1 , and define 



fn 



n 

n — L 



< n < \ - 1 
k <n< L— 1 



(3.3) 



The wavenumber corresponding to index «, is k\ = 2%f n jL. 
Transverse-smoothing is achieved by making an exponential 
suppression of the large wavenumber modes in the directions 
/ , m which are transverse to the direction j, and transforming 
the result back to position space. The prescription is to modify 
the Fourier-transformed array elements by the replacement 



1 



exp 



p(4+/i)k(ni,n 2j tt 3 ) (3-4) 



where p is a variational parameter, followed by an inverse 
Fourier transform back to position space. Then 



Aj(x,t,j) = fl(«l,« 2 ,«3) 



(3.5) 



is the field variable A? on a time-slice f, with suppression 
of the large wavenumber components in the /^-directions 
transverse to j. From this, we define the transverse-smoothed 
matrix-valued field variable 



Ai(x,t,j) = A?{x,t,j) 



In the same way, define 
Bi(x,t) 



l-±Tr[[/ ; (x,0 



(3.6) 



(3.7) 



Then, just as with the AJ'(x,f), suppress high-wavenumber 
components in directions transverse to some direction j, as 
described above. Denote the resulting field as B,(x,f,y). It is 
also convenient to define for ; ^ j, the average 



Bi(x,tJ) 



UAi{x,t,j)+Ai(x-ei,t,j) 



(Bi(x,t,j)+Bi(x-e h t,j) (3.8) 



derived from two links in the i direction which touch the line, 
in the j direction, running from quark to antiquark. 



AUxj) 



(3.1) 



For each lattice we loop through time-slices t, space compo- 
nents i = 1,2,3, color components c = 1,2,3, and an addi- 
tional direction j = 1,2,3, which will be associated with the 
direction of a line through the quark and antiquark. At each 
i,j,c,t, and denoting x = («i,«2,"3), (with indices running 
from to L — 1 ) define the three-dimensional array 



fl(ni,«2,«3) = A;(x,f) 



(3.2) 



1 In the two-gluon states described below we have enlarged this interval 
slightly, allowing for x coordinates up to two lattice spacings outside the 
interval. 



Having obtained the transverse-smoothed field variables 
Ai(x,t,j), B[(x,t,j), the next step is to define the six oper- 
ators {Qk} from which we obtain six states 



l*> = e*|o) (£=i-6) 



(3.9) 



which are used to build the trial gluon-chain state. Let the 
antiquark and quark charges lie at lattice sites xo and xr — 
xq + Rej, respectively, where ey is a unit vector in the j- 
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direction. 2 Then we choose 
Gl(f) = 12 

R-l 

Q 2 (t) = £A y (x + ney,t,7) 

n=0 

23(0 = E L ^;'(xo + ney,f,;)Ay(xo + n'ey,/,;) 

n=— 2 n'=n 
tf+2 fi+2 

24(0 = £ £ EA^xo + He/^^KxoWej,?,;) 

n=— 2 n'=n i^j 
R-l 

Q 5 {t) = Y< B j(xo + 11*1,1,1)12 

n=0 
R-l 

26(0 = ££5 ; (xo + «e,,f,;)l 2 (3.10) 

Qi is the zero constituent gluon operator, Q2 is the one gluon 
operator (one power of A), and the (?3-6 are two-gluon opera- 
tors, containing two powers of the A-field. The corresponding 
set of states {\k)} are not orthogonal, but this is taken care of 
by the Gram-Schmidt orthogonalization procedure described 
in the last section. 

The above choice of operators {Qk} is dictated by simplic- 
ity (only one variational parameter), and a certain amount of 
trial and error. It is certainly possible to invent operators cre- 
ating more general and sophisticated trial states, at the cost of 
additional variational parameters. 



IV. RESULTS 

We have carried out the calculation outlined in section 2, 
with the operators Qk described in section 3, for SU(2) lat- 
tice gauge theory at coupling p = 2.2 on a 12 4 lattice volume, 
P = 2.3 on a 16 4 lattice volume, and p = 2.4 on a 22 4 lattice 
volume. The operators Qk introduced in the previous section 
contain a variational parameter p, and we have calculated ma- 
trix elements 7y and V c hain(R) f° r eacn R, at twelve values of 
P 

p„ = (n-l)Ap (l<n<12) (4.1) 

with Ap = 0.025 at p = 2.2, and Ap = 0.02 at p = 2.3,2.4. 
The choice of n which minimizes V c hain(R) of course depends 
on both p and the quark separation. For example, the best 
choice at p = 2.2 and quark separation R = 1 was « = 2. At 
P = 2.4 and R = 9, the optimal value was n = 8. The results 
reported below are those obtained at the optimal value of p in 
each case. 

One matrix element of T which does not depend on the 
variational parameter is the zero-gluon to zero-gluon matrix 



2 Of course, in evaluating the transfer matrix elements T mn by lattice Monte 
Carlo simulation, we average over all xq and j = 1,2,3. 



element 

7li = <<>i|T|<t>i> = <1|^|1) (4.2) 

The color Coulomb potential reported in refs. [1, 2] is derived 
from this matrix element, i.e. 

V r COM /(/?) = -log(7i 1 ) (4.3) 

For purposes of comparison, we have also computed the usual 
static quark potential V true (R), which is the minimal possible 
field energy of a static quark-antiquark system, by the stan- 
dard method of computing timelike Wilson loops with "fat" 
spacelike links, and looking for a plateau in the lattice loga- 
rithmic time derivative, as described in, e.g., ref. [7]. 



A. Potentials 

The results for the potential V c i u ,i n (R) extracted from the 
variational state, compared to the Coulomb potential V C0M /(J?), 
and the static quark potential V trU e {R) computed by standard 
methods, are shown in Fig. 2. In this and all subsequent fig- 
ures, statistical errors are at most as large as the symbol sizes. 
The first point to note is that inclusion of one and two gluon 
operators alters the slope, but not the linearity, of the poten- 
tial. The second point is that the additional one and two-gluon 
operators in the trial state gives an estimate for the asymp- 
totic string tension which is much closer to the true value. We 
extract string tensions, in each case, from a fit to 

v(R)=cR-^+c (4.4) 

and find that while the Coulomb string tension differs from 
the true string tension by a factor of 2.3 (at p = 2.4), the string 
tension of V c hain{R) differs by roughly 38%. The lower figure 
is still a significant discrepancy, but it is also a considerable 
improvement over the zero-gluon Coulomb result. 

The one and two-gluon trial-state operators, described in 
the previous section, are not the only ones that can be imag- 
ined, and it may be that some modest improvement in these 
operators would bring string tension derived from the varia- 
tional state much closer to the true value. 



B. Gluon content 

The energy expectation value of the one-gluon state, ex- 
tracted from the one-gluon to one-gluon element of the trans- 
fer matrix, is estimated from the one-gluon to one-gluon ele- 
ment of the transfer matrix 

(/?) = - log r 22 (4.5) 

At small quark separation, this energy greatly exceeds the 
Coulomb potential V cou i, so the one constituent gluon content 
of the minimal energy state is expected to be much smaller 
than the zero constituent gluon content. However, as R in- 
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FIG. 3: Energy vs. quark separation of the zero-gluon and one-gluon 
states at (3 = 2.4. 



creases, this situation changes. In Fig. 3. we compare the 
energies of the zero and one-gluon states (|(j)i) and optimal 
\§2) states, respectively), as a function of quark separation R, 
at P = 2.4. In this case R is given in physical units, using 
the usual conversion from lattice spacing to fermis with string 
tension o = (440 Mev) 2 . We see that the Coulomb energy 
(energy of the zero gluon state) rises to meet V\ (R) at roughly 
one fermi. It is reasonable to expect that the one-gluon content 
of the minimal energy will rise accordingly. 
Let 



|V(*)> = £«*(*)H>*) 

k=\ 



(4.6) 



(c) 



denote the eigenstate of largest eigenvalue X max of the 6x6 
transfer matrix Ttj ,in the basis This is the minimal 

energy variational state. By zero and one-gluon "content" of 
this state, we mean the squared overlap | (<)),• |\|/)| 2 for i = 1,2 
respectively. Thus a\ gives the fraction of the norm due to the 
zero-gluon state, a 2 the fraction due to the one-gluon state, 
and 1 — a 2 — «2 is the fraction due to the two-gluon states, 
orthogonal to (j)i.2- We have plotted these fractions vs. R in 
physical units, for data at p = 2.2,2.3,2.4, with the result is 
displayed in Fig. 4. The gluon content vs. R in physical units 
is almost coupling independent. This scaling is gratifying, and 
a serves as a check of the whole procedure. It should be noted 
that the one-gluon content of the minimal energy state rises to 
equal the zero-gluon content at about one fermi, which is also 
the distance where the energies of the zero and and one-gluon 
states are roughly equal. 



FIG. 2: The color Coulomb potential V coll i(R), the "gluon-chain" 
potential V c i, a j n (R) derived from the variational state, and the static 
quark potential Vi„ie(R) extracted from "fat-link" Wilson loops. Re- 
sults are shown at lattice couplings (a) (3 = 2.2; (b) (3 = 2.3; and (c) 
(3 = 2.4. Continuous lines are from a fit of data points to eq. (4.4). 



C. Coulomb energy and finite size effects 

As mentioned in the Introduction, a strong objection to 
effective one-gluon exchange models of confinement is that 
such models lead inevitably to long-range dipole interactions. 
We may see a hint of this "dipole problem" in examining the 
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FIG. 4: Zero, one, and two-gluon content (fraction of the norm of the 
variational state) vs. quark separation R in fermis, at (3 = 2.2, 2.3, 2.4. 



lattice volume dependence of the Coulomb potential. In Fig. 
5 we plot V co „i(R) and V c i u ,i„(R) vs. R in lattice units, for the 
coupling (3 = 2.4 and lattice volumes 12 4 ,16 4 ,22 4 . On all 
three lattice L 4 lattice volumes, we see that the color Coulomb 
potential seems to bend away from linear at R m L/2. This de- 
parture from linearity is clearly a finite size effect, and is only 
seen in the vicinity of the largest possible on-axis quark sepa- 
rations. 3 However, it is very interesting to compare the sensi- 
tivity of the color Coulomb potential to finite-size effects, with 
the absence of such effects in the data for V c h a in(R)> which is 
also shown in Fig. 5. 

The different sensitivities of the Coulomb and chain po- 
tentials to finite lattice size could well be associated with the 
dipole problem. In particular, the color Coulomb field asso- 
ciated with the zero-gluon quark-antiquark state |(j)i) is not 
expected to be collimated in a flux tube. If the color-electric 
Coulomb dipole field extends throughout the lattice volume 
for the largest quark separations, then it is reasonable to ex- 
pect a corresponding sensitivity to the finite lattice volume. A 
cutoff in lattice volume would cut out a non-negligible part 
of the long-range field, and hence a portion of the Coulomb 
energy. By contrast, if the color-electric field of the minimal 
energy state |\|/) is largely collimated in a region whose trans- 
verse dimensions are small compared to the lattice size, then 
even for R w L/2 one would expect much less sensitivity to 
the finite lattice volume. In fact, in Fig. 5, we see no finite- 
size sensitivity whatever in V c i M j„(R), and the chain potential 
remains linear out to the largest on-axis separations at each 
lattice volume. This suggests that the chain state has no long- 
range dipole field, or at least that the long-range field is greatly 
suppressed relative to the color dipole field of the zero-gluon 
state. 



3 Our previous work on the color Coulomb potential [ 1 ] made use of off-axis 
separations for R ps L/2, which alleviated the distortion due to finite size 
effects. 



lattice volume dependence, (5=2.4 




2 4 6 8 10 12 

R 

FIG. 5: Sensitivity of the Coulomb potential V cou i(R) (solid sym- 
bols), and insensitivity of the chain potential V c f, a i n (R) (open sym- 
bols), to lattice volume. Data is for the gauge coupling P = 2.4, and 
lattice volumes L 4 = 12 4 , 16 4 ,22 4 . Quark-antiquark separation R is 
in lattice units. 



V. CONCLUSIONS 

The color Coulomb potential is the energy of a zero 
constituent-gluon state in Coulomb gauge, with two static 
quark-antiquark charges. The potential of such a state is 
known to rise linearly, at large quark separations, but its slope 
is far higher than the known asymptotic string tension. What 
we have found in this article is that the inclusion of a few 
gluon constituents, in simple trial states, does not alter the lin- 
ear rise of energy with quark separation, but does bring that 
energy very much closer to that of the true static quark poten- 
tial. According to our data, at one fermi we are only beginning 
to see the formation of a gluon chain state. At that distance the 
energies of the zero and one-constituent gluon states are about 
equal, and the zero and one-gluon states contribute equally to 
the minimal energy state containing a static quark-antiquark 
pair. 

Many questions remain. First, to what extent is the long- 
range dipole problem eliminated by inclusion of a few con- 
stituent gluons? We have seen an indication that the dipole 
problem is greatly reduced in the gluon-chain state, as com- 
pared to the zero-gluon state, but this needs further confir- 
mation. Secondly, since there is nothing compelling about 
the choice of operators in section 3, is it possible to con- 
struct a different set of operators and corresponding basis 
states such that the energy of the gluon-chain state 

comes significantly closer to the static quark potential? Fi- 
nally, it would be interesting to investigate chain formation 
for quark-antiquark sources in higher representations, and in 
larger gauge groups. For example, for N-ality k = 2 sources, 
do the quarks pull out a single chain of gluons, or two sep- 
arate chains, as they move apart? Are the results consistent 
with, e.g., Casimir scaling, or some other rule? We hope to 
return to these issues in a later publication. 
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